Historic,  Archive  Document 

Do  not  assume  content  reflects  current 
scientific  knowledge,  policies,  or  practices. 

- 


Numerical  Simulation  of  Snow  Avalanche 

T.  E.  Lang  and  R.  I.  Brawn 


Abstract 


A  computer  model  to  predict  snow  avalanche  impact  on  rigid  wall 
structures  is  developed  based  on  hydrodynamic  equations  for 
viscous  flow,  and  admitting  a  frictional  slip-plane  lower  boundary. 
Slope  normal  and  vertical  wall  configurations  which  extend  higher 
than  the  height  of  avalanche  flow  are  evaluated  for  normal  and 
shear  force  and  overturning  moment  from  avalanche  impact.  Impact 
forces  and  moment  time  depend  upon  avalanche  leading  edge  shape. 
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Numerical  Simulation  of  Snow  Avalanche 

Impact  on  Structures  c,  ^ 
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^ang  and  R.  L.^Brown 


Management  Implications 


Power  transmission  line  towers  and  bases,  avalanche 
sheds,  flow  diversion  barriers,  and  energy  dissipative 
structures  of  various  shapes  are  types  of  structures 
subject  to  snow  avalanche  flow  and  impact.  Design  of 
this  type  of  structure  has  developed  largely  from 
experience  and  includes  large  factors  of  safety,  which 
boost  material  and  labor  costs.  Little  engineering  infor- 
mation is  available  to  aid  in  establishing  expected 
forces  and  stresses,  particularly  in  situations  involving 
snow  impact. 

This  paper  presents  results  from  two  analyses  of 
structures  impacted  by  avalanches  at  various  speeds. 
Forces  on  the  structures  are  presented  in  graphical 
form,  so  that  design  specialists  may  refer  to  these  data 
to  estimate  forces  and  pressures  for  structures  that 
correspond  to  those  treated  in  the  report. 
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The/computer  programjdeveloped  to  obtain  the  data 
described  above  is  a  general-purpose  program  avail- 
able for  use  by  persons  involved  in  structural  design. 
Structures  of  many  different  shapes  may  be  input  and 
subjected  to  snow  avalanches  of  different  speeds  and 
depths.  Output  from  the  code  are  contact  pressures 
and  forces,  flow  profiles,  and  flow  velocities  of  the 
snow  as  it  impacts  and  flows  around  the  structure. 
Thus,  an  impact  becomes  a  pictorial  display  of  the 
processes  going  on. 

The  numerical  approach  to  representing  avalanche 
flow  and  impact  described  in  this  report  is  a  new 
development  in  the  methods  of  dealing  with  the  complex 
problems  of  snow  avalanche  mechanics.  Use  of  the 
program  by  others  is  encouraged  so  that  updates  and 
changes  can  be  made  based  upon  user  needs. 


Introduction 

When  snow  avalanches  impact  structures,  large 
forces  and  pressures  are  set  up  that  warrant  consider- 
ation in  structural  design.  Only  limited  data  exist  on 
the  forces  and  pressures  that  can  be  expected  from 
avalanches  of  different  speeds  and  depths  of  flow. 
Data  are  from  either  actual  physical  observations  of 
post-avalanche  impacted  structures,  or  from  planned 
experiments  over  limited  ranges  of  the  primary  param- 
eters of  speed-at-impact  and  flow  depth. 

Basic  contributions  to  knowledge  of  snow  impact 
were  made  by  Voellmy  (1964),  who  evaluated  evidence 
from  avalanches  in  Austria.  He  concluded  that  pres- 
sures in  the  range  of  10  to  20  t/m2  were  necessary  to 
account  for  the  damage  noted,  and  that  this  pressure 
range  is  attributable  to  the  solid  core  of  flowing  snow, 
and  not  to  the  airborne  component.  For  calculation  of 
impact  pressure,  Voellmy  proposed  a  version  of  the 
dynamic  pressure  equation  with  pressure  proportional 
to  the  square  of  the  flow  velocity  (or  kinetic  energy). 

Furukawa  (1957)  reported  results  of  releasing  snow 
blocks  to  slide  down  a  slope  and  impact  a  wall.  Average 
pressure  of  impact  was  measured  and  correlated  to 
other  parameters  by  the  equation 
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applicable  in  the  range  of  densities  150<e<550  kg/m3 
and  impact  velocity  5<v<20  m/s.  The  time  constant  of 
the  recording  oscillograph  of  these  experiments  is 
0.01  s,  so  that  rapid  variations  in  the  loading  are 
smooth.  Pressures  computed  by  Furukawa's  equation 
fall  into  the  same  range  as  reported  by  Voellmy. 

Saito  et  al.  (1963)  reported  the  results  of  experiments 
made  on  various  avalanche  control  structures.  Ava- 
lanche impact  forces  on  posts  set  up  in  the  avalanche 
path  were  measured.  Maximum  pressure  was  30  t/m2 
with  a  range  downward  to  an  average  of  around  20  t/m2. 
No  information  is  given  on  depth  of  flow,  nor  on  the 
nominal  density  of  the  impacting  snow. 

Salm  (1964)  presented  results  similar  to  those  of 
Furukawa,  except  that  force-recording  equipment  had 
a  shorter  time  constant,  so  that  millisecond  duration 
pulses  were  detected.  For  snow  blocks  impacting  at 
velocities  of  11  to  13  m/s,  average  pressures  were 
obtained  in  the  range  15  to  20  t/m2  for  impact  onto  an 
elastically  soft  wall.  Maximum  pressures  of  load  pulses 
of  several  milliseconds  duration  were  recorded  and 
found  to  be  two  to  five  times  larger  than  average  pres- 
sures reported  above. 

Perhaps  the  most  complete  set  of  data  on  avalanche 
impact  is  reported  by  Schaerer  (1973)  based  upon  a 
number  of  avalanche  measurements  made  at  Rodgers 
Pass,  British  Columbia,  Canada. 
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Pressure  gages  mounted  on  posts  in  avalanche  paths 
were  used  to  monitor  average  impact  pressure.  Meas- 
urements also  were  made  of  depth  of  flow,  nominal 
density  of  the  snow  debris  and  average  avalanche 
speed.  Based  upon  a  dynamic  pressure  equation, 
Schaerer  compared  computed  and  measured  pressures. 
Average  pressures  ranged  from  3  to  26  t/m2  for  velocity 
variations  from  15  to  53  m/s  and  flow  depths  between 
1.5  and  1.8  m  as  evidenced  from  snow  deposits  on  trees 
and  side  walls  of  the  flow  channels. 

Additional  evidence  is  reported  by  Mears  (1975) 
based  upon  analysis  of  damaged  trees  in  the  path  of  an 
avalanche  in  Colorado.  For  flow  speeds  in  the  range  of 
18  m/s,  loading  was  estimated  at  8  to  10  t/m2  for  a  flow 
height  of  1.1  m  and  snow  density  of  300  kg/m3. 

All  of  the  above  results  tend  to  indicate  average 
dynamic  pressures  up  to  30  t/m2  for  flow  velocities  up 
to  40  m/s;  although  incomplete  data  reporting  tends  to 
prevent  conclusive  ranging.  Shimuzu  et  al.  (1973) 
reported  different  results  based  upon  measurements  of 
three  avalanches  in  Kurobe  Canyon,  Japan.  Indications 
from  the  discussion  in  this  report  are  that  the  ava- 
lanches were  "high  speed;"  however,  no  information 
was  presented  on  the  actual  range  of  velocities  of 
the  avalanches.  Pressures  were  reported  from  32  to 
134  t/m2  depending  upon  type  of  recording  system,  and 
a  maximum  pressure  of  210  t/m2  is  mentioned  without 
elaboration  or  explanation.  Lacking  information  on 
recording  system  response  and  velocities  of  the  flows, 
it  is  possible  to  attribute  these  high  pressures  to  either 
high  impact  velocity,  or  to  reporting  of  maximum  pres- 
sures of  load  pulses  rather  than  of  pressures  averaged 
over  10  to  100  ms  or  more,  as  is  done  by  Salm  (1964) 
and  Schaerer  (1973). 

Having  an  indication  of  the  nature  of  the  force  (dis- 
tribution in  time)  of  typical  avalanche  impact,  the 
possibility  of  computer  simulation  of  the  phenomenon 
was  investigated. 

This  report  summarizes  new  methodology  for  pre- 
dicting impact  of  snow  avalanches  on  structures.  The 
approach  involves  computer  modeling  of  the  flow  and 
impact  using  an  iterative  solution  of  the  Navier-Stokes 
equations  of  motion  for  viscous,  Newtonian  fluid  flow. 
Time  and  spatial  advance  of  the  fluid  through  the  finite 
difference  grid  of  the  flow  domain  is  carried  out 
by  iterative  refinement  of  the  momentum  and  mass 
distributions  consistent  with  the  governing  equations. 
Original  reporting  of  a  computer  algorithm  for  two- 
dimensional  flow  in  either  a  closed  or  free-surface 
domain  is  by  Hirt  et  al.  (1975).  Modification  of  this 
general  purpose  code  to  the  specific  modeling  of 
avalanche  flow  is  reported  by  Lang  et  al.  (1979),  which 
includes  considerable  detail  on  the  characteristics  and 
application  of  the  code  to  avalanche  flow  problems. 
This  report  further  modifies  this  code  to  facilitate 
calculation  of  pressure  and  force  distributions  on  rigid 
structures  subjected  to  avalanche  impact.  Computer 
simulation  data  are  compared  with  existing  physical 
information  to  evaluate  its  accuracy.  Physical  proper- 
ties assumed  for  snow  in  the  computer  model  are 
nominal  properties  that  have  been  determined  to  be 


representative  for  avalanche  runout  prediction  (Lang 
et  al.  1979). 

Computer  simulation  could  be  used  to  develop  design 
criteria  for  structures  without  having  to  accumulate 
physical  data  from  random  avalanches. 

Governing  Equations  and  Flow  Parameters 

The  flow  characteristics  of  a  possible  airborne  or 
snow-dust  component  of  an  avalanche  must  be  distin- 
guished from  the  flow  of  the  denser  or  core  component. 
Our  consideration  is  directed  to  the  core  component, 
which  has  a  typical  velocity  range  of  20  to  50  m/s. 
Whether  or  not  the  core  material  behaves  the  same 
when  surrounded  by  an  airborne  component  is  un- 
known. For  snow  having  large  viscosity,  the  flow  is 
assumed  to  be  incompressible,  which  is  a  conservative 
approximation  when  considering  pressures  and  forces 
of  impact  upon  obstacles. 

The  numerical  analysis  of  snow  avalanche  impact  is 
carried  out  for  a  "typical"  cross-section  of  the  flow,  so 
that  it  is  sufficient  to  consider  two-dimensional  forms 
of  the  governing  equations.  For  an  incompressible, 
viscous  fluid,  the  equations  of  motion  are  the  Navier- 
Stokes  equations: 

du  _  _  1_  _dg_  f  d*u  d2u] 
dt  q    dx  +v\dx2  +  dy2) 

[1] 

dv  _y_]_  _9p_  f  d2v  d2v] 
dt  q    dy  +P[dx2  +  dy2  J 

where 

u,v  =  components  of  velocity  (m/s) 

t  =  time  (s) 
X,Y  =  components  of  of  body  force  per  unit 
of  density  (m2/s2) 
q  =  fluid  density  kg/m3 
p  =  fluid  pressure  (Pa) 
v  =  kinematic  viscosity  (m2/s) 
x,y  =  rectangular  Cartesian  coordinates  (m) 

and 

d        d  d  d 

dt  ~  dt  +U  dx  +V  dy 

Equation  [1]  is  not  sufficient  to  solve  for  the  pressure, 
p,  and  two  components  of  velocity,  u  and  v.  The  addi- 
tional equation  needed  is  the  mass  continuity  equation, 
which  in  two  dimensions  for  either  homogeneous  or 
heterogenous  incompressible  flow  takes  the  form 


This  completes  the  number  of  equations  needed  to 
solve  for  the  unknowns  of  the  problem. 
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The  stresses  acting  upon  a  fluid  element  are: 

Tyy=-p  +  2M-|^  [3] 
dv  du 

where  \i  is  the  dynamic  viscosity,  which  is  related  to  the 
kinematic  viscosity  by  v  =  \aIq.  The  pressure,  p,  is 
assumed  to  be  compressional,  and  hence  a  negative 
contribution  to  stress.  Equation  [3]  is  needed  to  com- 
pute the  forces  that  develop  on  the  impacted  obstacle 
as  a  result  of  fluid  contact. 

Finite  difference  versions  of  equations  [1]  and  [2]  are 
actually  used  in  computer  code  AVALNCH  (Lang  et  al. 
1979)  to  simulate  avalanche  flow  and  impact.  The  code 
is  written  so  that  variable  kinematic  viscosity  and 
surface  friction  coefficients  can  be  used  to  represent 
the  flow.  Flow  down  a  slope  is  simulated  by  assigning 
gravity  components  to  the  body  force  intensities,  X  and 
Y,  of  equation  [1].  Actual  computed  parameter  values 
from  the  code  are  normalized  by  material  density,  so 
that  computed  pressure  is  actually  in  units  pig  (m2/s2), 
stress  is  t^q  (m2/s2),  and  total  impact  force  is  FJq 
(m4/s2).  The  finite  difference  grid  that  represents  the 
flow  domain  plus  a  one-cell  boundary  layer  is  shown  in 
figure  1. 


Figure  1.— Arrangement  of  finite  difference  grid,  avalanche 
leading  edge,  and  slope-normal  impact  wall. 


Figure  2.— Boundary  cell  representation. 


Angle  4>  is  the  angle  of  the  slope  on  which  the  ava- 
lanche advances.  Angle  if>L  is  the  initial  angle  of  the 
leading  edge  of  the  avalanche,  which  is  assumed  to  be 
straight.  Parameters  6X  and  6y  are  the  x  and  y  dimen- 
sions, respectively,  of  an  individual  cell  of  the  uniform 
grid.  Surface  friction  is  defined  by  specifying  the  x 
component  of  velocity  in  the  lower  boundary  cell  using 
the  Newtonian  shear  law.  If  u2  is  the  x  component  of 
velocity  in  the  lowest  flow  domain  cell,  and  ui  is  the 
corresponding  velocity  in  the  lower  boundary  cell  (fig. 
2),  then 

U,=U:  (1-2/)  [4] 

where  /  is  the  friction  coefficient. 
If  /  =  0,  then  Ui  =  u2  and  the  boundary  is  slip-free. 
If/ =  1,  then  iii  =  -  Uz  and  velocity  at  the  interface  be- 
tween the  cells  is  zero,  which  is  a  no-slip  condition.  The 
friction  coefficient  is  related  to  shear  stress  at  the 
interface,  as  follows: 

Finally,  friction  coefficient,  /,  is  assigned  to  each  lower 
boundary  cell  individually,  so  that  friction  may  vary 
along  the  slope,  depending  upon  physical  circumstances 
of  the  problem. 

In  the  remainder  of  this  paper,  set  v  =  0.5  m2/s  and 
/  =  0.5,  as  representative  of  nominal  values  for  these 
parameters. 
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Impact  Upon  a  Slope-Normal  Wall 


Problem  Definition 

Much  of  the  experimental  data  from  load  cell  meas- 
urements are  taken  from  test  structures  in  which  the 
sensing  elements  are  arrayed  in  a  line  normal  to  the 
slope  surface.  Thus,  direct  comparison  of  experimental 
data  and  computer  simulation  data  is  possible  using  a 
slope-normal  wall  (i.e.,  a  wall  perpendicular  to  the 
local  slope).  This  geometry  is  also  the  simplest  case  to 
model  using  computer  code  AVALNCH,  since  the  wall 
representation  amounts  to  a  zeroing  of  velocity  com- 
ponents in  a  single  column  of  cells  (fig.  1). 

Computer  definition  of  the  problem  is  shown  in  figure 
1.  The  avalanche  leading  edge  is  initially  defined  by 
angle  \pL,  and  material  in  all  cells  to  the  left  of  the 
leading  edge  are  given  an  initial  slope-parallel  speed  in 
the  range  of  equilibrium  speed  for  avalanche  flow  on  a 
slope  of  angle  <}>.  This  initial  assumption  is  not  correct; 
but  in  the  first  cycle  of  computations  to  establish  mass 
continuity,  the  velocity  and  pressure  values  in  each  cell 
are  adjusted  to  meet  specified  compatibility  require- 
ments. This  results  in  a  change  of  shape  of  the  leading 
edge  profile  from  a  straight  line  to  other  shapes, 
examples  of  which  are  given  later  as  part  of  the  pres- 
entation of  data.  If  the  initial  assumption  of  slope- 
parallel  speed  is  close  to  the  nominal  speed  for  a  slope 
of  angle  <J>,  and  with  the  assumed  values  for  surface 
friction  and  viscosity,  then  the  compatibility  correction 
is  small.  What  ensues  is  a  flow  in  which  the  leading 
edge  changes  shape,  but  for  which  the  nominal  speed 
of  the  entire  mass  does  not  systematically  change  with 
advance  of  the  flow.  In  all  cases,  computer  runs  were 
made  with  different  initial  velocity  estimates  until  the 
flow  showed  the  persistence  in  velocity  as  described 
above. 

Under  the  driving  attraction  of  gravity,  the  ava- 
lanche flows  down  the  slope,  impacts  the  obstacle,  and 
piles  up  against  it.  The  numerical  algorithm  used  in 
program  AVALNCH  does  not  admit  representation  of 
flow  over  the  top  of  the  obstacle,  so,  consideration  is 
limited  to  initial  buildup  of  the  impacting  snow.  A  more 
intricate  computer  code  that  admits  the  overflow 
condition  is  being  adapted  to  avalanche  flow. 


Force  Computations 


Figure  3.— Parameters  and  geometry  of  cells  at  the  impact  wall. 
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The  upper  limit  on  the  sum,  n,  is  the  number  of  cells  in 
the  itix  column  that  contains  snow.  For  the  top  cell  that 
contains  snow  it  may  be  necessary  to  use  a  partial 
height  rather  than  6y  in  equation  6. 

For  cell  (i,  j)  the  normal  force  on  the  wall,  using  equa- 
tion 3,  is: 


An  important  consideration  is  the  definition  of  the 
resultant  forces  that  act  on  the  impacted  wall,  as  com- 
puted from  velocity  and  pressure  data  determined  for 
each  cell.  A  free  body  diagram  of  the  normal  and  shear 
forces  that  act  between  a  cell  in  the  itil  column  and  jth 
row  and  the  wall  is  shown  in  figure  3.  Also  defined  are 
the  cell  velocity  components  and  the  resultant  base 
forces  FN  and  Fg,  and  base  moment,  Mg,  needed  for 
equilibrium  at  the  wall.  In  the  case  of  the  slope-normal 
wall,  for  the  ith  column  of  cells, 


where  (Pjj/e)  is  the  actual  quantity  computed  as  "pres- 
sure" in  AVALNCH.  The  forces  reported  are  also  scaled 
to  unit  density  by  equation  [7]  so  that  the  units  of  the 
term  in  equation  [7]  are  (m4/s2),  which  assumes  force  is 
per  unit  distance  in  the  third  coordinate  direction. 
Total  normal  force  on  the  wall  is  then 
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1  =  2 


Various  options  can  be  exercised  in  computing  shear 
stress  (and  also  shear  force).  One  possibility,  which  is 
based  on  the  assumption  of  a  no-slip  wall,  and  a  local 
finite  difference  definition  of  the  derivatives  for  shear 
stress  in  equation  [3]  is: 
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Recognizing  zero  terms  and  cancelations  in  this  equa- 
tion, it  reduces  to 
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=  y  (Uj _  u  +  1-ui_lj_1}- v{vLj  +  Vj  j _  x)  [9] 


Alternately,  by  a  coarser  grid  averaging  in  which  a 
single  shear  stress  is  defined  for  each  cell,  the  following 
shear  force  for  cell      is  defined 


L  =  "4  (Uj-ij  +  i-Ui-ij-i)  -Jj(vi-itj  +  v1_1(M)  [10] 
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Although  a  number  of  finite-difference  methods  can  be 
used  to  define  shear  stress,  equation  [9]  is  used  in 
subsequent  evaluations.  Using  equation  [10]  implies  a 
surface  friction,  /,  of  0.5,  whereas  equation  [9],  with 
the  no-slip  condition,  implies  /=  1.0.  Total  shear  force 
on  the  wall  is  then 


I 

1  =  2 


Q 


[11] 


Summary  and  Evaluation  of  Data 
for  a  Slope-Normal  Wall 

Setup  of  the  problem  is  shown  in  figure  1.  The  total 
flow  domain  is  modeled  by  30  cells  along  the  slope,  and 
30  cells  normal  to  the  slope.  All  cells  are  squares  0.05  m 
on  a  side.  The  slope-normal  wall  is  placed  in  cell  number 
24,  and  the  leading  edge  of  the  avalanche,  at  the  start 
of  flow,  is  in  cell  number  21. 

Output  from  program  AVALNCH  consists  of  a  listing 
of  the  pressure  per  unit  density,  velocity  components, 
and  flow  height  for  each  cell  at  any  specified  time 
during  the  impact.  Forces,  FnIq  and  FsIq,  moment,  MbIq, 
and  maximum  pressure,  pmIq,  are  computed  at  each 
time  increment  of  the  computations  and  are  listed 
sequentially.  A  typical  impact  condition  is  considered 
to  outline  the  results.  Slope  angle  was  set  at  4>  =  30°, 
leading  edge  angle  i^L  =  30°,  nominal  slope-parallel 
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Figure  4.  — Distribution  of  snow  upslope  from  a  rigid  slope-normal 
wall  of  several  times  measured  from  time  of  impact  at  T  =  0.0  s. 
Slope  angle  =  30°.  Impact  velocity  =  20  m/s.  Kinematic 
viscosity  =  0.5  m2/s.  Friction  coefficient  =  0.5. 

flow  velocity  U0  =  20  m/s,  and  viscosity  and  friction  at 
v  =  /  =  0.5.  One  type  of  result  is  the  change  of  shape  of 
the  leading  edge  of  the  avalanche  as  the  flow  pro- 
gresses (fig.  4).  At  the  instant  of  impact,  the  leading 
elements  of  the  flow  have  changed  from  an  initial  angle 
of  ^»L  =  30°,  to  an  impact  angle  if»j  =  76°.  After  initial 
contact,  the  avalanche  front  piles  against  the  slope- 
normal  wall  with  a  profile  sequence  as  depicted  in 
figure  4.  Maximum  pressure  and  maximum  normal 
force  on  the  wall  occur  when  the  flow  surface  is 
approximately  parallel  to  the  slope  (at  T  =  8.3  ms,  for 
this  case).  The  time  histories  of  the  normal  and  shear 
forces  and  the  moment  at  the  base  of  the  wall  are 
shown  in  figure  5.  It  is  seen  that  a  sharp  peak  in  normal 
force  occurs  at  approximately  8.3  ms  into  the  impact, 
and  persists  for  approximately  1.5  ms. 

The  distribution  of  pressure  along  the  slope-normal 
wall  for  the  various  times  into  the  impact,  corres- 
ponding to  the  different  snow  profiles  of  figure  4,  are 
shown  in  figure  6.  The  pressure  distribution  corres- 
ponding to  T  =  8.3  ms  is  0.1  ms  off  the  time  at  which 
the  pressure  is  maximum  for  the  entire  impact,  at 
pM/@  =  2,500  m2/s2.  Maximum  pressure  was  found  to 
vary  with  the  angle  of  the  leading  edge  of  the  ava- 
lanche, ipL.  As  ipL  is  changed,  the  impact  angle  ipj  also 
changes  according  to  the  results  listed  in  table  1. 
Although  leading  edge  angle  is  more  than  doubled, 
impact  angle  changes  by  only  13°,  but  with  the  angle 
approaching  90°.  Thus,  the  slope  of  the  frontal  face  of 
the  avalanche  at  impact  has  a  pronounced  affect  on 
maximum  pressure  and  total  force  exerted  on  the  wall. 
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Figure  5.— Force  and  moment  components  as  functions  of  time 
after  impact  upon  a  slope-normal  wall.  Slope  angle  =  30°. 
Initial  leading  edge  angle  =  30°.  Impact  velocity  =  20  m/s. 
Kinematic  viscosity  =  0.5  m2/s.  Friction  coefficient  =  0.5. 


Q. 

CD 
O 

c 
o 

Q 


0.80 

0.70 
0.60 

0.50 
0.40 
0.30 
0.20 
0.10 

0 


T=8.3ms 

)  ^? 

i 

i               i                i  i 

500         1,000         1,500  2,000 
Pressure/unit  density  (m2/s2 ) 


2,500 


0.02  0.04  0.06 

Time  after  impact  (s) 


Figure  7.— Pressure  profile  at  height  of  maximum  pressure  on 
wall  versus  time  after  impact  for  avalanche  leading  edge  angles 

of  30°,  45°,  60°,  and  75°. 

Using  the  largest  value  of  maximum  pressure  per 
unit  density  of  9,700  m2/s2  (table  1),  convert  this  to  an 
actual  pressure  value  by  assuming  a  nominal  snow 
density  of  300  kg/m3  so  that 


=  ePMH^oo-^ 


=  2,910,000  kg/ms2  or  N/m2 
=  (2,910,000  N/m2)  (1  kg/9.807N) 
(1  ton/1,000  kg) 
297  ton/m2 


Figure  6.— Distribution  of  pressure  per  unit  density  along  slope- 
normal  wall  at  several  times  measured  from  time  of  impact  at 
T  =  0.0  s. 


Pact  ~ 


which  is  approximately  a  factor  of  10  greater  than 
experimentally  measured  pressures  which  are  in  the 
range  of  30  t/m2  as  noted  in  the  introduction. 

From  the  physical  measurements  reported  by  Schea- 
rer  (1973),  maximum  pressure  is  associated  with  time 
durations  on  the  order  of  0.1  second.  This  occurs 
because  of  the  relatively  long  time  constant  of  the 
load-cell  type  sensor  used  to  monitor  force.  Thus  a 
pressure  peak  of  1.5  ms  duration  would  not  be  detected 
except  in  an  averaging  sense.  Therefore,  it  is  necessary 
to  consider  the  pressure  distribution  over  0.1  second 
duration  in  the  analytical  results  in  order  to  make 
comparison  with  experimental  data.  To  this  end,  the 
pressure  distribution  at  the  height  of  maximum  pres- 
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Table  1.— Variation  in  maximum  pressure  and  forces  as  a  function  of  leading  edge  angle,  xpL, 

for  the  slope-normal  wall  geometry 


Maximum  pressure         Maximum  normal        Maximum  shear       Maximum  bending 


Leading  edge 

Estimated  impact 

per  unit 

force  per 

force  per 

moment  per 

angle,  ipL 

angle,  t^, 

density 

unit  density 

unit  density 

unit  density 

degrees 

degrees 

m2ls2 

m4ls2 

/777S2 

m5/s2 

30 

76 

2,500 

830 

200 

240 

45 

81 

4,700 

1,750 

320 

500 

en 
OU 

81 

6,700 

2,800 

500 

1  ,uuu 

75 

86 

8,700 

5,800 

1,300 

3,600 

80 

88 

9,700 

7,200 

1,600 

5,000 

Table  2.— Average  pressure  on  the  wall  at  the  height  of  maximum  pressure  for  different  leading 

edge  angles  and  for  a  nominal  flow  velocity  of  20  m/s 

Leading  edge 

Average  pressure 

Average 

Height  on  wall 

angle 

per  unit  density1 

pressure2 

of  maximum  pressure 

Maximum  pressure2 

degrees 

m2ls2 

tonlm2 

m 

tonlm2 

30 

305 

9.3 

0.53 

76.5 

45 

455 

13.9 

0.55 

144 

60 

600 

18.4 

0.63 

205 

75 

815 

25.0 

0.95 

266 

80 

850 

26.0 

1.25 

297 

'Based  on  0.1 -s  time  duration  since  impact. 

2Based  upon  snow  density  of  300  kglm3. 

sure  on  the  wall  for  four  different  cases  is  shown  in 
figure  7.  From  these  data  the  average  pressure  can  be 
determined  from  the  area  under  each  curve.  The  com- 
puted values,  and  the  conversion  to  actual  pressures, 
assuming  a  snow  density  of  300  kg/m3,  are  listed  in 
table  2. 

By  averaging  pressure  over  0.1  s,  it  is  seen  that 
computed  pressures  fall  in  the  same  range  as  experi- 
mentally measured  pressures.  Plotting  maximum  and 
average  pressures  as  functions  of  the  estimated  impact 
angles  (fig.  8],  linear  variation,  particularly  of  average 
pressure,  is  a  good  approximation.  It  should  be  noted 
that  computation  of  impact  angle  is  approximate  in 
that  the  0.05-  by  0.05-m  grid  assumed  in  determining 
these  results  is  not  adequate  for  detailed  resolution  of 
the  leading  edge  shape.  Because  of  this  limitation, 
impact  angle  is  based  upon  straight  line  approxima- 
tions of  leading  edge  shape,  as  depicted  in  figures  4, 
A-l,  A-3,  and  A-5. 

By  changing  the  angle  of  the  ground  surface,  ava- 
lanches with  different  nominal  flow  velocities  can  be 
modeled.  By  doing  this,  flow  velocities  of  U0  =  7,  10,  15, 
and  30  m/s  were  modeled  at  initial  leading  edge  angles 
of  30°,  45°,  60°,  and  75°.  For  any  given  initial  leading 
edge  angle,  the  impact  angle  for  all  flow  speeds  was 
found  to  be  approximately  the  same  as  that  reported 
for  the  20  m/s  case  (table  1).  Thus,  using  different  values 


12,000- 


10,000 


CM 

cn 


~  8,000 


-g  6,000  - 


£  4,000 

CO 
CD 


2,000- 


80  82  84 
Impact  angle  (degrees) 


Figure  8.— Variation  in  maximum  and  average  pressures  with 
estimated  leading  edge  angle  for  imapct  at  20  m/s  onto  a  slope- 
normal  wall.  Kinematic  viscosity  =  0.5  m2/s.  Friction  coeffi- 
cient =  0.5. 
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of  leading  edge  angle  at  impact,  the  variation  in  maxi- 
mum and  average  pressure  versus  the  nominal  flow 
speed  of  the  avalanche  is  shown  in  figure  9.  Results 
obtained  for  U0  =  7  m/s  show  no  pressure  peak,  but 
rather  a  monotonic  increase  in  pressure  as  the  snow 
piles  against  the  wall. 

A  surface  undulation  of  the  snow  profile  after  impact 
subsides  approximately  0.2  m  upslope  from  the  im- 
pacted wall  (fig.  4).  Farther  upslope  from  this  station, 
the  particles  at  the  surface  advance  at  the  nominal 
speed  of  20  m/s.  It  can  be  hypothesized  that  the  height 
of  snow  at  the  terminal  point  of  the  undulation  is  repre- 
sentative of  the  height  needed  to  produce  the  maximum 
(or  average)  pressure  of  the  impact,  dependent  also 
upon  impact  angle.  These  observed  flow  heights  are 
summarized  for  the  different  nominal  flow  speeds  and 
impact  angles  in  figure  10.  Combining  the  results  of 
figures  9  and  10  led  to  the  conclusion,  for  example,  that 
an  avalanche  moving  at  20  m/s  and  having  a  nominal 
leading  edge  angle  at  impact  of  i^j  =  76°,  will  produce 
an  average  pressure  of  300  m2/s2  per  unit  of  density,  a 
maximum  pressure  of  approximately  2,500  m2/s2  per 
unit  of  density,  and  an  equivalent  flow  height  of  0.55  m. 
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Figure  9.— Variation  in  pressure  versus  nominal  flow  velocity  for 
different  leading  edge  angles  at  impact  with  a  slope-normal 
wall. 


Figure  10.— Equivalent  flow  height  versus  nominal  flow  velocity 
for  different  leading  edge  angles  at  impact  with  a  slope-normal 
wall. 


Impact  Upon  a  Vertical  Wall 
Problem  Description 

The  vertical  wall  case  represents  the  geometry  when 
an  avalanche  impacts  against  a  building  or  similar  type 
structure.  Again,  only  a  two-dimensional  approximation 
of  the  flow,  on  what  would  likely  be  a  three-dimensional 
problem  unless  the  structure  and  avalanche  are 
broad,  is  considered. 

As  in  the  slope-normal  case,  the  avalanche  is  de- 
picted by  a  constant  slope  leading  edge  starting  in  cell 
number  21  of  a  30  by  30  grid.  The  case  of  a  nominal 
flow  velocity  of  20  m/s  is  considered  to  compare  max- 
imum pressure  and  average  pressure  with  the  slope- 
normal  case.  If  differences  are  not  large,  then  the  cost 
of  duplicate  computer  studies  of  the  vertical  wall  case 
are  not  necessary,  and  reference  can  be  made  to  the 
results  of  the  section  on  governing  equations  and  flow 
parameters  as  typical  for  both  cases.  A  nominal  flow 
velocity  of  20  m/s  corresponds  to  equilibrium  flow  on  a 
slope  of  30°.  Thus,  to  specify  a  vertical  wall,  zero 
velocity  cells  must  be  specified  in  an  array  as  shown  in 
figure  11.  Starting  at  the  bottom  surface  of  the  grid  the 
2-2-1  sequencing  of  steps  toward  the  uphill  averages  to 
an  equivalent  vertical  wall. 

Listing  of  the  computer  code  for  the  vertical  wall 
case  is  given  in  Appendix  B.  Substitute  instructions  are 
also  inserted  for  the  case  of  the  slope-normal  wall.  For 
problem  arrays  other  than  30  by  30  the  DIMENSION 
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statement  of  lines  9  and  10  must  be  changed.  For  a 
nominal  impact  speed  other  than  20  m/s,  instruction 
112  must  be  changed.  If  the  leading  edge  of  the  ava- 
lanche is  not  to  start  in  cell  21,  then  instruction  110 
changes.  All  remaining  data  for  this  code,  particularly 
the  input  data,  is  the  same  as  described  for  program 
AVALNCH  (Lang  et  al.  1979).  Specifically,  input 
instructions  are  13,  19,  65,  66,  and  68,  for  which 
description  of  the  parameters  and  arrays  are  given  in 
Appendix  A  (Lang  et  al.  1979). 

Force  Calculations 

For  the  vertical  wall  geometry  it  is  necessary  to 
distinguish  between  an  open  cell  and  a  corner  cell 
(fig.  11).  For  the  open  cell,  the  calculation  of  normal  and 
shear  force  at  the  interface  with  the  wall  is  that  speci- 
fied by  equations  7  and  9.  For  a  corner  cell,  two  normal 
and  two  shear  forces  are  specified  on  the  two  adjacent 
fixed  boundaries  (fig.  12).  If  the  cell  number  is  i,j  cor- 
responding to  the  column  and  row  number,  respectively, 
then  the  normal  forces  are  computed  from  the  equations 

where  the  subscripts  follow  the  order  as  depicted  in 
figure  3.  The  two  shear  forces  for  the  cell  are  com- 
puted by  the  equations: 


Figure  11.— Grid  array  and  wall  geometry  for  vertical  wall  impact. 
Case  shown  is  a  30°  slope,  for  which  nominal  flow  speed  is 
20  m/s  and  leading  edge  angle  is  ipL  =  30°. 


i 


Figure  12.— Forces  acting  on  the  boundary  of  corner  cell  i,  j. 


Fs2 

—  =  -  yfai-ij  +  Vi-ij-i) 

where  6X  =  6y  has  resulted  in  simplification  of  these 
equations. 

By  summing  the  respective  slope-parallel  and  slope- 
normal  forces  of  each  cell  over  all  cells  that  contain 
snow,  the  total  slope-parallel  and  slope-normal  forces 
are  obtained.  Total  moment  at  the  base  of  the  wall  must 
necessarily  include  the  moment  of  all  force  components 
Fjv  ,  Fsj,  Fjv2,  and  F$2  of  a  corner  cell,  and  Fjv  and  F$ 
of  an  open  cell  as  the  cells  are  offset  in  moving  vertically 
upward. 

Summary  and  Evaluation  of  Data 
for  a  Vertical  Wall 

The  pressure  variations  at  the  height  of  maximum 
pressure  for  the  four  cases  of  leading  edge  angles  of 
xpL  =  30°,  45°,  60°,  and  75°  are  shown  in  figure  13.  The 
evaluation  is  for  the  case  of  a  nominal  impact  velocity 
Uo  =  20  m/s,  kinematic  viscosity  v  =  0.5  m2/s,  friction 
at  /  =  0.5,  and  a  slope  angle  <j>  =  30°.  These  values  are 
selected  so  that  the  flow  in  the  mainstream  of  the  ad- 
vancing avalanche  is  equilibrium  flow. 

For  the  case  of  the  vertical  wall  geometry,  the  angle 
between  the  wall  and  the  upslope  ground  surface  is 
60°.  For  the  avalanche  cases  considered,  leading  edge 
angles  increase  rapidly,  and  all  angles  at  the  instant  of 
impact  for  the  slope-normal  geometry  equal  or  exceed 
\pj  =  76°.  At  these  angles  of  the  leading  edge  in  the  case 
of  a  vertical  wall,  upper  reaches  of  the  leading  edge 
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will  contact  the  wall  before  the  base.  This  is  in  contrast 
to  the  occurrence  in  the  slope-normal  wall  geometry  in 
which  initial  contact  is  at  the  base.  Thus,  different 
pressure  and  force  histories  can  be  expected  for  the 
two  geometries.  Maximum  pressures  and  average 
pressures  over  the  initial  100-ms  duration  of  impact  as 
functions  of  leading  edge  angle  are  given  in  table  3  for 
the  vertical  wall  geometry.  Although  peak  pressures 
show  a  wide  variation  among  the  different  cases, 
average  pressure  varies  by  less  than  a  factor  of  2.0. 
This  is  verified  by  the  near  equality  of  areas  under  the 
pressure-time  curves  of  figure  13. 

Detailed  variations  in  normal  force,  shear  force,  and 
bending  moment  on  the  wall  for  the  case  of  the  30°, 
leading  edge  angle  avalanche  are  shown  in  figure  14. 


Similar  plots  for  the  cases  of  leading  edge  angles  of 
45°,  60°,  and  75°  are  presented  in  Appendix  A  (figs. 
A-7,  A-8  and  A-9).  No  attempt  is  made  at  plotting  the 
snow  surface  distribution  for  the  vertical  wall  case 
because  of  uncertainty  as  to  the  area  of  contact  at 
the  leading  edge  as  a  function  of  time.  To  accurately 
evaluate  contact  area  as  a  function  of  time  requires 
only  refinement  in  the  finite-difference  grid,  which  is 
technologically  straightforward.  Results,  at  least  for 
the  case  of  \pL  =  30°,  show  that  maximum  pressure  and 
force  on  the  vertical  wall  occur  when  the  impacting 
snow  front  piles  up  until  the  snow  surface  is  approxi- 
mately parallel  to  the  ground  surface  next  to  the  impact 
surface.  This  is  consistent  with  the  results  obtained 
also  for  the  slope-normal  wall. 


Table  3.— Variation  in  pressures  as  functions  of  leading  edge  angle  for  the  vertical  wall 

geometry 


Leading  edge 
angle 

Maximum  pressure 
per  unit  density 

Average  pressure 
per  unit  density1 

Average  pressure2 

degrees 

m2ls2 

m2/s2 

ton/m2 

30 

730 

280 

8.6 

45 

1,200 

320 

9.8 

60 

2,200 

340 

10.4 

75 

10,600 

450 

13.8 

1  Average  pressure  over  initial  100  ms  of  impact. 
2Based  upon  snow  density  weight  of  3,000  kg/m3. 
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Figure  13.  — Pressure  profile  on  vertical  wall  versus  time  after 
impact  for  avalanche  leading  edge  angles  of  30°,  45°,  60°,  and 
75°.  Nominal  flow  velocity  is  U0  =  20  m/s.  Kinematic  viscosity 
=  0.5  m2/s.  Friction  coefficient  =  0.5. 


Figure  14.  — Force  and  moment  components  as  functions  of  time 
after  impact  on  a  vertical  wall.  Slope  angle  =  30°.  Initial 
leading  edge  angle  =  30°.  Impact  speed  =  20  m/s.  Kinematic 
viscosity  =  0.5  m2/s.  Friction  coefficient  =  0.5. 
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Table  4.— Average  and  maximum  pressures  generated  on  slope-normal  and  vertical  walls  as  a 
function  of  leading  edge  impact  angle,  ipv  and  for  a  nominal  flow  velocity  of  20  m/s  at  impact 


 Slope  normal  wall    Vertical  wall  

Leading  edge  Maximum  pressure  Average  pressure1  Maximum  pressure  Average  pressure1 
impact  angle         per  unit  density  per  unit  density  per  unit  density  per  unit  density 


degrees  m2/s2  m2/s2  m2/s2  m2/s2 

76  2,500  305                               730  280 

80  4,700  455  1,200  320 

81  6,700  600  2,200  340 
86  8,700  815  10,600  450 
88  9,700  850                              —  — 


1  Averaged  over  initial  100  ms  of  impact. 


Conclusions  from  Slope-Normal  and 
Vertical  Wall  Results 

Although  the  computer  solutions  include  calculations 
of  forces  and  moments  on  the  slope-normal  and  vertical 
walls  (which  is  essential  information  for  design),  pres- 
sure is  the  fundamental  quantity  for  comparison  of  the 
two  cases.  The  average  and  maximum  pressures  per 
unit  density  for  the  two  cases  are  given  in  table  4. 

For  impact  angles  of  76°,  80°,  and  81°,  the  slope- 
normal  wall  pressures  are  larger  than  corresponding 
vertical  wall  pressures  with  larger  differences  noted 
for  maximum  pressures  than  for  average  pressures. 
For  an  impact  angle  of  86°,  the  vertical  wall  maximum 
pressure  suddenly  jumps  to  a  value  larger  than  the 
corresponding  value  for  the  slope-normal  case.  A 
plausible  explanation  for  this  is  that  the  leading  edge 
of  the  initial  impacting  snow  is  near-parallel  to  the  wall 
surface  at  contact.  The  result  of  this  is  a  generated 
maximum  pressure  of  10,600  m2/s2  per  unit  density  that 
peaks  and  subsides  in  less  than  a  millisecond  (fig.  A-9). 
This  is  followed  by  a  pressure  tailoff  such  that  average 
pressure  greater  than  100  ms  is  comparable  to  that  of 
the  other  impact  angle  cases.  This  wave  "slapping" 
action,  well  recognized  in  water  wave  phenomenon, 
apparently  can  occur  in  avalanche  impact  also,  and 
constitutes  the  extreme  case  in  pressure  generation. 

To  further  verify  this  occurrence,  a  fifth  case  of 
impact  on  the  slope-normal  wall  was  run.  An  initial 
leading  edge  angle  of  80°  changes  into  an  estimated 
impact  angle  of  88°  which  results  in  a  maximum 
pressure  of  9,700  m2/s2.  Assuming  a  straight-line 
extrapolation  on  maximum  (and  average)  pressure  ver- 
sus impact  angle  (fig.  8),  the  estimated  limit  pressure 
curves  of  figure  9  are  determined  for  the  case  of  a 
slope-normal  wall.  Although  not  explicitly  verified  for 
different  nominal  flow  speeds,  computed  results  indi- 
cate that  limit  pressures  for  the  vertical  wall  case  are 
of  the  same  order  as  for  the  slope-normal  case.  Since 
all  pressures  reported  include  static  pressure  plus  a 
velocity  gradient  term  (equation  [3]),  it  is  expected  that 
pressures  on  walls  that  make  small  angles  with  the 
ground  surface  would  be  considerably  smaller  than  the 
limit  pressures  reported  in  figure  9. 

It  is  noted  also  that  the  limit  pressures  cited  are 
based  upon  avalanche  impact  with  a  rigid  wall.  Any 


wall  flexibility  can  be  expected  to  decrease  the  maxi- 
mum pressure.  Additionally,  compressibility  of  the 
snow  itself  will  also  tend  to  decrease  pressure.  Neither 
wall  flexibility  nor  snow  compressibility  are  modeled  in 
the  computer  representation.  Note,  if  wave  reflection 
inside  a  specific  impacted  structural  system  is  con- 
sidered, internal  stresses  may  develop  that  are  larger 
than  the  incident  impact  pressure. 

A  more  systematic  development  of  pressure  peaking 
is  noted  in  the  slope-normal  case.  There  is  a  sharp  peak 
for  all  impact  angles  evaluated.  The  nominal  duration 
of  each  peak  is  from  1.0  to  2.0  ms.  Assuming  the 
existence  of  this  type  of  loading,  the  question  arises 
whether  a  pressure  pulse  of  this  duration  can  have 
significant  effect  on  structural  materials.  Results  of 
tests  on  small  steel  specimens  (Newmark  and  Halti- 
wanger  1962)  show  that  stresses  on  the  order  of  twice 
the  yield  stress,  then  applied  rapidly,  can  cause 
material  yielding  and  fracture  within  1.0  ms.  However, 
extrapolation  of  this  to  reinforcing  steel  is  not  appar- 
ent. Although  stresses  on  structures  may  develop  to 
values  beyond  yield  or  fracture  values,  strains  neces- 
sary for  the  failure  will  not  develop  because  of  large 
inertia  effects.  Thus,  for  typical  structures  made  of 
steel,  concrete,  or  wood,  it  is  unlikely  that  the  millisecond 
duration  pulse  will  have  a  significant  effect,  and  design 
should  be  predicated  upon  average  pressure  response. 

From  these  results  we  conclude  that  the  specific 
shape  of  the  leading  edge  of  the  avalanche  at  impact 
has  a  significant  effect  on  the  maximum  pressure 
that  is  generated.  A  worst  case  geometry  has  been 
evaluated  from  this  computer  study;  however,  the 
statistical  probability  of  its  occurrence  is  not  known. 
To  obtain  supportive  data  from  field  observations, 
advanced  avalanche  observation  techniques  would 
have  to  be  developed  that  would  allow  penetration  of 
the  usual  snow-air  cloud  that  accompanies  all  but  the 
slow-moving  wet-snow  avalanches. 

A  second  recourse,  although  less  accurate  than 
physical  measurements,  would  center  on  determining 
the  viscous  and  friction  parameters  of  the  flow  more 
accurately  than  what  is  currently  the  case.  Then,  using 
this  information,  steady  state  flow  would  be  simulated 
on  the  computer,  and  a  statistical  estimate  of  the 
distribution  of  leading  edge  angle  would  be  determined 
by  numerical  means. 
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A  third  (experimental)  alternative  is  to  use  a  number 
of  sensitive  pressure  transducers  with  millisecond 
response  capability  on  a  rigid  post  and  measure  the 
relative  times  of  pressure  buildup  between  transduc- 
ers. This  provides  a  means  for  estimating  the  leading 
edge  shape  of  the  more  dense  flowing  snow  of  an 
avalanche.  From  a  number  of  such  measurements,  a 
statistical  description  of  leading  edge  shape  would  be 
obtained. 

Until  such  time  that  more  data  is  available  on  ava- 
lanche flow,  recourse  to  worst  case  design  is  the  safest 
approach,  although  probably  not  very  economical. 

Finally,  the  maximum  pressure  peaks  reported  by 
Salm  (1964)  and  the  data  given  here  are  compared.  The 
peaks  measured  by  Salm  were  two  to  five  times  greater 
than  average  pressures  reported,  and  the  pulses  show 
durations  on  the  order  of  5  to  10  ms.  Results  reported 
here  show  maximum  pressure  to  be  as  much  as  nine 
times  greater  than  average  pressure  and  the  pulses  of 
duration  1  to  2  ms.  Whereas,  the  computed  results  are 
for  impact  upon  a  rigid  wall,  the  impact  in  Salm's 
experiments  is  against  an  elastic  wall.  Based  upon 
physical  reasoning,  one  would  expect  a  decrease  in 
pulse  height  and  spreading  of  the  pulse  when  impact  is 
changed  from  a  rigid  to  an  elastic  wall.  Qualitative 
verification  of  this  is  possible  in  program  AVALNCH 
by  inserting  a  cell  at  the  wall  over  which  velocity  is 
greatly  reduced  rather  than  set  to  zero.  Although  this 
has  not  been  done,  differences  in  initial  peak-load 
characteristics  between  the  computer  results  and 
those  of  Salm  seem  to  have  a  plausible  physical  basis. 

The  general  conclusion  from  these  computer  results 
is  that  program  AVALNCH  has  modeling  capability  for 
avalanche  impact.  Uncertainties  as  to  leading  edge 
geometry  prevent  explicit  ranging  of  expected  pres- 
sures; however,  preliminary  results  show  that  upper 
limits  on  pressure  can  be  established.  Further  cor- 
relation between  computer  predictions  and  field 
observations  are  warranted;  the  types  of  parameters 
that  need  to  be  measured  to  provide  closer  correlations 
have  been  identified. 
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Appendix  A 
Supplemental  Plotted  Data 
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Figure  A-1 .  — Distribution  of  snow  upslope  from  rigid  impact  wall 
at  several  times  measured  from  timeofimpact  at  T  =  0.0. 
Slope  angle  =  30°.  Leading  edge  angle  =  45°.  Flow  velocity  = 
20  m/s.  Kinematic  viscosity  =  0.5  m2/s.  Friction  coefficient 
=  0.5. 
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Figure  A-2.— Force  and  moment  components  as  functions  of  time 
after  impact  on  a  slope-normal  wall.  Leading  edge  angle=45°. 
Flow  velocity=20  m/s.  Kinematic  viscosity=0.5  m2/s.  Friction 
coefficient=0.5. 


13 


1.20  - 


Distance  from  wall  (m) 


Figure  A-3.— Distribution  of  snow  upslope  from  rigid  impact  wall 
at  several  times  measured  from  time-of-impact  at  T  =  0.0. 
Slope  angle  =  30°.  Leading  edge  angle  =  60°.  Flow  velocity 
=  20  m/s.  Kinematic  viscosity  =  0.5  m2/s.  Friction  coefficient 
=  0.5. 
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Figure  A-4.— Force  and  moment  components  as  functions  of  time 
after  impact  on  a  slope-normal  wall.  Leading  edge  angle  = 
60°.  Flow  velocity  =  20  m/s.  Kinematic  viscosity  =  0.5  m2/s. 
Friction  coefficient  =  0.5. 
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Figure  A-5.— Distribution  of  snow  upslope  from  rigid  impact  wall 
at  several  times  measured  from  time-of-impact  at  T  =  0.0. 
Slope  angle  =  30°.  Leading  edge  angle  =  75°.  Flow  velocity  = 
20  m/s.  Kinematic  viscosity  =  0.5  m2/s.  Friction  coefficient  = 
0.5. 
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Figure  A-6.— Force  and  moment  components  as  functions  of  time 
after  impact  on  a  slope-normal  wall.  Leading  edge  angle  =  75°. 
Flow  velocity  =  20  m/s.  Kinematic  viscosity  =  0.5  m2/s. 
Friction  coefficient  =  0.5. 
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Figure  A-7.— Force  and  moment  components  as  functions  of  time  Figure  A-8.— Force  and  moment  components  as  functions  of  time 

after  impact  on  a  slope-vertical  wall.  Slope  angle  =  30°.  after  impact  on  a  slope-vertical  wall.  Slope  angle  =  30°. 

Leading  edge  angle  =  45°.  Flow  velocity  =  20  m/s.  Kinematic  Leading  edge  angle  =  60°.  Flow  velocity  =  20  m/s.  Kinematic 

viscosity  =  0.5  m2/s.  Friction  coefficient  =  0.5.  viscosity  =  0.5  m2/s.  Friction  coefficient  =  0.5. 
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Figure  A-9.— Force  and  moment  components  as  functions  of  time 
after  impact  on  a  slope-vertical  wall.  Slope  angle  =  30°. 
Leading  edge  angle  =  75°.  Flow  velocity  =  20  m/s.  Kinematic 
viscosity  =  0.5  m2/s.  Friction  coefficient  =  0.5. 
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Appendix  B 


Listing  of  Computer  Code  for  Vertical 
Wall  and  Inserted  Slope-Normal  Wall 


C     PROGRAM  IMPACT  -  VERTICAL  WALL.  ^  [      SLOPE  NORMAL  WALL  j 

C       EQUIVALENT  HORIZONTAL  GRID  OPTION 

C  THIS  PROGRAM  HAS  BEEN  DEVELQPEO  TO  NUMERICALLY   SGLVE  (USING 

C  FINTTE  DIFFERENCE  TECHNIQUES)  THE  NAVTFR  -  STOKES  EQUATIONS 

C  FOR  TRANSIENT   FLUID  FLOW  PROBLEMS.     IT  WILL   BE  USEO  TC  SIMULATE 

C  THE  PROBLEM  OF  SNOW  AND  ICE   AVALANCHES  ON   SLOPES  OF   VARYING  SLOPE. 

DIMENSION  UC3Q,30),V(3Q,30),UN(30f  30).Vn(30.  30).PO0  ,30), 
1XPUTC  8),H(   30)fHN(  30), JT(  30), NAME(2C) t GX(   3C),GY(   30),FRC(  30) 

REAL  NU 

INTEGER  CYCLE 

READ  (105,45)  NAME 

WRITE  (108,35) 

WRITE  (108,45)  NAME 

C  *  *  READ  AND  PRINT  INITIAL   INPUT  DATA 
READ(105,25)  (XPUT( I) , 1=1 t 8 ) 

IBAR=XPUT(1)   ;  J8AR=XPUT(2)   ;  DELX=XPUT(3)   \  DELY=XPUT(4) 
NU=XPuT(5)     I     FRK=XPUT(6)     %     TWFIN=XPUT(7 )      I     C  WP  R  T  =XP  UT  (  8  ) 
WRITE  (108,50)  (XPUT(T),T=1,8) 

25  FORMAT(8F10.O) 

35  FORM AT( 1H1) 

45  FORMAT(20A4> 

47  FQRMAT(6X*I*7X*J'8X'U'13X'V'13X'P'13X'H'9Xf'SUR  CELL*) 

48  FORMAT(4X1T3,5X,T3,4(4XilPE10.3),6Xf 12) 

49  F0RMAT(1X'CYCL='I3, IX 'ITR= '13 f 1 X *TME = * 1PE9 .2 , IX *SNH= 'E 9 .2 , 1 X ,  'UM=# 
1E9.2,1X*FX=*E9.2, IX 'F Y= *E 9. 2 , 1 X *T B= *E 9 . 2 , 1 X * P M^ *E 9. 2  ) 

50  FORM ATC1H   ,1X*IBAR= 'F4. 0. 2X *JB AR=*F3.0 ,2X'DELX='F6.2,2X'DELY='F5.2 
1.2X'NU='F4.2,2X'FRK='F4.2,2X*TWFTN='F5.0  9Z  X  'C  WPRT  =  'F5  .  0  ) 

60  FORMAT(SFIO.O) 

61  FORMAK8F10.3) 

70  F0RMAT(1H0,35X.'FL0W  HEIGHT*) 

71  FOR  MAT( 1H0 , 25X  EL  EV  AT  TON  CHANGE  FOR  EACH  CELL*) 

72  FORMAT(lH0,25X, 'BOUNDARY  FRICTION  COEFFICIENTS*) 

73  FqRMAT(1H0,25X, 'SLOPE-PARALLEL  GRAVITY  COMPONENT*) 

74  F0RMAT(1HQ,25X, 'SLOPE-NORMAL  GRAVITY  COMPONENT') 

75  FORMAT(lH0,30X,*END  OF  INPUT  DATA*//) 

82  F0RMAT(5X, 'PROBLEM  RUNNING  TIME  EXCEEDED-C ALCUL AT  IONS   TERMINATED  *) 

* 

C  *  *  COMPUTE   CONSTANT  TERMS  AND  INITIALIZE   NECESSARY  VARIABLES 


IMAX=IBAR+2   ;  JMAX=JBAR+2 

TM1=TMAX-1      I  JM1=JMAX-1 

RDX=1.3/DELX 

RDY=1.0/DELY 

DELM  =  CE  LY/100 • 

DELT=.0005 

IM2=lMAX-2 

JM2=JMAX-2 

T=FLG=Q.O 

CYCLE=ITER=INO=0 

G=9.8  ;   OMG=1.7  ;   EpSI=0.5      ;   ALPHA=0.1    ;    GAMMA=C.l    ?  DZR0=1.Q 

BET  A=  0MG/(2.*DELT*(R0X**2+RDY**2  )) 

TCPRT=INT(CWPRT) 

IF(ICPRT.EQ.I)  ICPRT=2 

DO  100  I=1,IMAX 

H(I)=HN(I)=JT(I)=GX(I)=GY(I)=FRC(  I)=0.0 
DO  100  J=1,JMAX 
100  U(I,J)=V(I, J)=UN(I, J)=VN(If J)=P(I, J)=0.a 

C  *  *   SPECIAL   INPUT  DATA 


READ(10  5,60)(H(T),T=2  ,TM1) 
READ(10  5,60XHN(I),I=2,IMI) 
IF(FRK.GT.O.O)  GO  TO  120 
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READU05.60)  (FRC(I )f 1=2, IM1) 

GO  TO  130 
120  DO  125  T=2,TM1 
125  FRC(I)=FRK 
130  CONTINUE 

DO  150  1=2, IM1 

SP=HN(I)/DELX 

CP=SQRT(1.0-SP*SP> 

GX(I*=G*SP 
150  GY(I)=-G*CP 

WRITE(1G8,70) 

WRTT£a08,61)(H(T^f  T=2f  THl) 
WRITE(108,7l5 

WRITE C  108,  61KHNCD  ,1=2  tIMl) 
WRITE(108,725 

WRITE(108,61)  CFRCCI),T=2,TM1) 
WRITEC108.73) 

WRITEU08,61KGX(n,I=2,IMl) 
WRITE(108,74> 

WRITEUG8, 61  KG Y(I), 1  =  2  ,1*1) 
WRTTE(108,75) 

C       DETERMINE   INITIAL  TOP   SURFACE  INDICES. 

DO  240   1=2, IM1 
JT(I)=INT(H(I)*R0Y+l.E-6)+2 
IFCjTCD.GT  .J»U)  JT(I)=JM1 
240  HN(I)=0.0 
H(l)=  H(2) 
H(IMAX)=  H(IMi) 
JTC1)=  JT(2> 
JT(I«AX)=  JT(IKl) 

* 

C  *  *  CALCULATE  HYDROSTATIC  PRESSURE 


DO  280  1=2, IM1 
jTl=JTCn 


DO  280  J=2,JT1 
280  PCT  ,  J)=-GY(I)*(H(  I)-(FLOAT(  J)-1.5  )*D£LO 

* 


C  *  *  SET   INITIAL  VELOCITY 

DO  310  T=2,21 
DO  310  J=2,JT(I) 
U(I,J)=2G. 
310  V(I,J)=0.0 

ASSIGN  4280   TC  KRET 
GO  TO  2000 

C  *  *  START  CYCLE 
* 

1000  CONTINUE 
TTER=0 
FLG=1. 

ASSIGN  3000  TO  KRET 

C  *  *  COMPUTE  TEMPORARY  U  AND  V 
♦ 

DO  1100  1=2, I"l 
JTl  =  JTCO 
DO  1100  J=2,JT1 


FUX=((UN( I, J) *UN( 1+1. UN(I9 J3+UNCI+1 ,J))+ ALPH A*A3S(UN(I  , J)+UN( 
1T+1 ,J))*CUN(Tf J)-UNCT^1, j))-(UN(T-l, J)+UN(T  f J))*(UNCT-1, J)+UN(I f J) 
2)-ALPHA*ABS(UNCI-l, J)+UN(I, J))*CUNCI-1, J)-UN(I,J)))/(4  .*DELX) 

FUY=((VN(I,  J>+VNCI*lf  J))*CUN(I,  J)  ♦UNCI,  J*l  )) 
14 ALPH A* ABS(VN(I.J)+VN(I+1 , J ) )* (UN ( I , J )-UN ( I f J  + 1 )  ) 

2-  CVNCT, J-l)+VN(f+l, J-1>)*CUN(T,J-1)+UN(T , J) 5 

3-  ALPHA*A*jSCVNCI,  J-1)+VN(I+1,  J-1))*(UN(  I ,  J- 1  )-UN(  I  ,  J  )  ) )  / (4  .  * 06  L Y  ) 
FVX=((UN(T, J)+UN( I, J+1))*(VN( I  , J)  +  VN(I  +  1 , J >)+ ALPH A* Ab S ( UN (  I  ,J)*UNC 

II,  J+1)>*(VN<I,  J)-VNCI*1  ,  J))-CUN(  1-1.  J)+UNCI-1,  J^1))*CVNCI-1,  J)^V;MC 
2I,J))-ALPHA*ABSCUN<I-1, J)+UN(I-lf  J  +  1))*(VN(I-1, J)-VN(I j J)>)/C4.*0E 
3LX) 

FVY=((VN(T, J)*VN(T, J*1))*(VN(T , J)+VN(T  §J  +  1))*ALPHA*A4S(VN(T ,  J)  + VN 
1(1, J*1))*(VN(I, J)-VN(I, J+1))-(VN( I, J-1)+VN( I , J))*( VNCI, J-l  >*VN(I,  J 
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1100 


-1)-VNCI» J)>>/<4.*DELY) 
'  ^/DELY**2 


,/DELX**2+ 
ut  j-1))/0ELY**2) 
1,J))*RQX   +  GXCIJ-FUX-^UY+VISX) 
J+1))*RDY   4  GY(I)-FVX-FVY+VTSY) 


C 
C 


*  *   SET   BOUNDARY  CONDITIONS 


2000  CONTINUE 

HNCl)=  HNC2) 
HNCX*AX)=HNCTM1) 
JTC1)=JTC2) 
JT£INAX)=JTCIM1) 
LEFT  VI ALL  CONTINUOUS  INFLOW. 
RIGHT  WALL  CONTINUOUS  OUTFLOW. 


IF(ITX 

DO  2250  J=1,J«AX 
UC1, J)=UC2, J) 
VC1. J)=VC2,J) 
UCTM1,J)=UCTM2,J) 
YCIMAX, J)=VCI*1,J) 
2200  CONTINUE 
C       TOP  WALL  CONTINUOUS 
C       BOTTOM  WALL  RIGID  - 
DO  2500  T=1.TNAX 
IFCITER.GT.O)  GO  TO  2400 
VCI,JMl)=VCI,JM2> 
UCI, JHAX)=UCI,JM1) 


R.GT.O)  GO  TO  2200 


OUTFLOW. 
WITH  FRICTION 


C 


2400 
2500 

♦  * 


* 

C 


vci,n=o.o 

UCI,1)=UCI, 2)*C1.0-2.0*FRCCI)) 

FREE   SURFACE  8 OUNDARY  CONDITIONS 

DO  2650   1  =  2, IM1 
JT1=JTCI) 

IFCJTCI+1).LT. JTCI))  UC I , JT 1 )  =  UC I , JT1- 1 ) 
VCI,  JT1)=  VCI, JT1-1 )-DELY*RDX*CUCI» JT1)- 
UCI,JT1U)=  UCItJTl) 


2650 

♦  *   SPECIAL  BOUNDARY  CONDITIONS 


UCI-1,JT1)) 


JNX  =  1 
J  =2 
1=24 

2900   IFCJNX.EQ.6)   JNX=1    ;   I  =  T-1 
IFCJNX.EQ.5)  1=1-1 
IFC JNX.EQ.3)  1=1-1 
IFCJ.EQ.21)  GO  TO  2920 
DO  2910  N=I,IM1 
UCN,J)=0.0 
2910  VCN+1,J)=0.0 
J=J*1 
JNX= JNX+1 
GO  TO  2900 
2920  CONTINUE 


 GO  TO  KRET, (3000,4280) 

3000  CONTINUE 


DO  2900  1=24, IM1 
00  2900  J=1,JH1 

UCIt J)=UO 
2900  V(I41,J)=0.Q 


* 

C  *  * 
* 


HAS  CONVERGENCE  BEEN  REACHED 

IFC FLG. EQ.C • )   GO  TO  4000 
ITER  =ITER+1 

IFCITER.LT. 500>  GO  TO  3050 
IFCCYCLE.LT. 10)  GO  TO  4000 
T=1.E+1Q 


GO  TO  4000 
3050  FLG=0.0 


JU^W      I   l_     —  J  •  u 

C  *  *  COMPUTE  UPDATED  CELL  PRESSURE  AND  VELOCITIES 
♦ 


JB1  =  2 
DO  3500   1=2, INI 
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JT1=JTCI> 
CO  3500  J=2,JT1 
TFCJ.EQ.JTl)  GO  TO  3100 
GO  TO  3200 

3100  ETA=0ELY/<HNCI)-(FL0ATCJT1)-2.5)*DELY) 
DELP  =  a.O-ETA)*PCl,  JT1-D-PCI,  JT1) 
GO  TO  3300 

3200  D=RDX*<UCI, J)-U(I-l,J))*RDr*<vCI, J)-V( I, J-l» 

IF(A8S(D/DZR0).GE.EPSI)  FLG=1.0 

DELP  =  -BE  T A*D 
3300  PCI,J)=PO, J)*DELP 

UCI, J)=U(I, J)+DELT*RDX*DELP 

UCI-1, J)=U(I-1,J)-DELT*RDX*D£LP 

VCI,J)=VCIjJ)+DELT*RDY*DELP 
3500  VCI, J-l)=V<It J-1)-DELT*RDY*0ELP 

GO  TO  2000 
4000  CONTINUE 

♦ 

C  *  *  COMPUTE   NEW   POSITION  FOR  TOP  SURFACE 

oo  4100  1=2, m 

JT1=JT(I) 

HV=   RDf *CHN(T )-FLOAT( JT1-2)*0ELY) 
UAV  =  0.5*(U(I-1,JT1)  ♦  UCItJTD) 

H(I)=HN(I)  +DELT*CHV*V(I»  JT1 )+( 1. O-HV )* V( I , J Tl 

1-1)-0.5*ROX«CUAV*HNCI+1)+GAMHA*ABSCUAV)*CHNCI)-HNCI+1 )) 
2       -UAV*HN(T-1)-GAMMA*ABS(UAV)*<HN<T-1)-HN(T )))) 

4100  CONTINUE 

* 

C  *  *  CALCULATE  CELL  IN  WHICH  SURFACE  IS  LOCATED  ANO  UPDATE  ARRAY 
* 

DO  4250   1=2, IM1 

JTCI)=INTCH(I)*RDY+1.0E-6)+2 

IFCjTCD.GT.  JMl)  JTCI)=JM1 
4250  CONTINUE 

ASSIGN  4280  TO  KRET 

GO  TO  2000 
4280  CONTINUE 

C  *  *  CALCULATE  TOTAL  FLUIO  VOLUME 
* 

FVOL=0.0 
00  4300   1=2, 1*1 
4300  FVOL=FVOL*HCI)*OELX 

* 

C   *  *   ADVANCE   U,V,H  ARRAYS. 

4910  UMAX=VMAX=O.Q 

DO  4900  I=1,INAX 

HNCI)=HCI) 

DO  4900  J=1,JMAX 

UNCI,J)=U(I,J) 
4900  VN(T  ,  J)  =  V(T , J) 

DO  4950   1=2, IM2 

DO  4950  J=2,JTCI) 

UT=ABS(UN(I,J)) 

VT=ABS(VN(I, JJ) 

TFCUT.GT.UMAX)  UMAX=UT 
4950   IF(VT.GT.VWAX)  VMAX=VT 
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THESE  INSTRUCTIONS 
REPLACED  BY  THOSE 
LISTED  ON  PAGE  21 . 


c  *  * 

_*  


COMPUTE  BASE  FORCES   AND  MOMENT  ON  WALL 


FX=Fr=Tb=PM=0.0 

JCY=1 

JNX  =  1 

j=j+i 

ifc jcy. eq.5 )  go  to  4990 

PA=PCIt J>+2.*NU*U(I-1, J^/OELY 
IFCPA.GT.PM)  PM=PA 

TFCPCT, J).GT.0.0. AND. JNX.EQ.l)  GO  TO  496* 


4955 


IFCPCl, J).GT.0.0. AND. JNX.EQ.2)  GO  TC  4965 
IFCPCl, J).GT.Q.Q. AND. JNX.EQ. 3)  GO  TO  496C 
IFCPCl, J).GT.5.0.AND. JNX.EQ.4)  GO  TO  496* 


IFCPCl, J). GT. 0.0. AND. JNX. EQ.5)  GO  TO  4965 
IFCJNX.EQ.l)  JNX  =  2 


IFCJNX.EQ.2)  JNX  =  3 
IFC  JNX. EQ. 3 )  JNX=4 
IFCJNX.EQ.4)  JNX=5 


GO  TO  4955 

1=1-1  ;   GO  TO  4955 

GO  TO  4955 

1=1-1  J   GO  TO  4955 

T=T-1  ;  JCr=Jcr+i  ; 


GO   TO  4955 


4965 


TFCJNX.EQ.5)  JNX=1 
4960  RX=PCI,J)*DELY*2. *NU*UC 1-1 , J) 
FX=FX*RX 

FY=FY*NU*CUCI-l,J*l)-UCI-l,J-l))/2.-NU*CvCI,  J)*V(I,  J-l)) 
TB=TB4RX*CFL0ATCJ)-1.5)*DELr 
JNX=JNX+1 
GO  TO  4955 

RX=PCI, J)*DELY*2.*NU*UCI-1,J> 
RY=PCI, J)*0ELY+2.«NU*VCI, J-l) 
SX=-NU*(V(I-1,  J-D  +  UCI-1,  J)  )/2. 
SY=-NU*CUCI-1, J-D  +  VCI, J-l))/2. 
FX=FX*RX*SX 
FY=FY*RY+SY 

TB=TB+RX*CFL0ATCJ)-1.5)*DELY+SX*CFL0ATCj)-l. ) *DE L r-R T * ( FLO A T ( T ) -2 
14.5)*DELX 
JNX=JNX*1 
1=1-1 

IFCJNX.EQ.6)  JNX=1   ;  JCY=JCY+1 
GO  TO  4955 
CONTINUE 

SNHI=H(24>  "  —  ■  


,4590 


C  *  * 
♦ 

5000 


5250 

♦ 

C  *  * 
6000 


LIST  VELOCITY,  PRESSURE,   AND  SURFACE  POSITION 

WRITE (108, 49)   CYCLE  .T T ER , T • SNHT , U * AX , FX , F Y , T B  , P M 

IF(CYCLE.EQ.ICPRT)  GO  TO  5030 

IFCCYCLE.NE.ICPRT)  GO  TO  6000 

ICPRT=  ICPRT  ♦  INTCCWPRT) 

CONTINUE 

WRTTE(l08,47) 

DO  5250  1=15,25 

JT1=  JTCI) 

JT2=JT1+1 

DO  5250  J=1.JT2 

WRTTEC108,48)  T , J ,U (T , J) , V(T , J ) , P ( I ,  J ) , H ( I )  ,  JT  1 
CONTINUE 

IFCIN0.EQ.2)  GO  TO  6520 

RECOMPUTE  CONTROL  PARAMETERS. 

IFCCYCLE.EQ.O)  GO  TO  6300 
DTX=DELX/UMAX 
DT Y=DELY/VM  AX 
DELT=AMIN1(DTX ,DTY)/3.Q 
IF(ITER.LT.IO)  DELT=1.5*DELT 
IFCNU-1.E-6.LT. 0.0)  GO  TC  6300 
OET=(DELX*DELY5**2/ (2 .*NU*(DELX**2+0ELY**2 )) 
IFCDELT.LT. DET)  GO  TO  6300 
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6300 


C  *  * 

6400 

6440 

6500 

6520 
6600 


DELT=0.9*DET 
T=T+DELT 

IF(CYCLE.EQ.O)  GO  TO  6400 
OAX=UKAX*DELT/DELX 
DAY=VHAX*DELT/DELY 
ALPHA=1.35*AMAX1(0AX,0AY) 
IFCALPH A.GT.l.O)  ALPHA=0.95 
G AHM  A= ALPHA 

8ETA  =  CMG/C2.*D£LT*CRDX**2+RDY**2  )) 
TEST   FOR   PROGRAM  TERMINATION, 

IFCT.GT.TWFIN)  INO=2 
IF(IND.GT.l)  GO  TO  6500 
CYCL£=CYCLE*1 
GO  TO  1000 
T=T-DELT 


GO  TO  5060 
WRITEC108, 
STOP 
END 


82) 


THESE  INSTRUCTIONS  FOR  THE  SLOPE-NORMAL  WALL  REPLACE  BOXED  INSTRUCTIONS  ON  PAGE  20. 


IF(H<24).GT.DELM)  GO  TO  4960 
FX=FY=TB=PH=0. 0 
GO  TO  4990 
4960  JU=INT<H<24)*R0Y*l.E-6)+2 
J=2 

FX=FY=TB=PH=0 .0 
4970  TF<J.EQ.JU)  G3  TO  4980 

RX=P<24,J)*DELY-2.0*NU*(U<24, J)-U(2  3, J)) 
FX  =  F  X—R  X 

FY=FY*NU*<U<23#J>D-UC2  3f  J-l))/2.0-NU*CV(24f  JJ  +  VC24,  J-U> 
TB=T8+RX«CFLOATCJ)-1.5)*OELY 

GO  TO  4970 
4980  HH=H(24)-DELY*FL0ATCJ-2) 

RX=CP<2  4, J)-2.*NU*<U(24f J)-U(23, J))/DELX)*HH 
PX=FX—RX 

FT=FT4NU*< (UC 23 t J )-UC 23, J- 1) ) /DE L Y-2.*V( 24, J- D/CELX) *HH 
TB=T84RX*CHC24)-HH/2. 0) 
00  4985   1=2. JU 

PA=PC24t T>*2.0*NU*U<2  3,T)/DELY 
4985   IFCPA.GT.PM>  PM=PA 
4990  CONTTNUF 
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U.S.  Department  of  Agriculture 
Forest  Service 

Rocky  Mountain  Forest  and 
Range  Experiment  Station 


The  Rocky  Mountain  Station  is  one  of  eight 
regional  experiment  stations,  plus  the  Forest 
Products  Laboratory  and  the  Washington  Office 
Staff,  that  make  up  the  Forest  Service  research 
organization. 

RESEARCH  FOCUS 

Research  programs  at  the  Rocky  Mountain 
Station  are  coordinated  with  area  universities  and 
with  other  institutions.  Many  studies  are 
conducted  on  a  cooperative  basis  to  accelerate 
solutions  to  problems  involving  range,  water, 
wildlife  and  fish  habitat,  human  and  community 
development,  timber,  recreation,  protection,  and 
multiresource  evaluation. 

RESEARCH  LOCATIONS 

Research  Work  Units  of  the  Rocky  Mountain 
Station  are  operated  in  cooperation  with 
universities  in  the  following  cities: 

Albuquerque,  New  Mexico 
Bottineau,  North  Dakota 
Flagstaff,  Arizona 
Fort  Collins,  Colorado* 
Laramie,  Wyoming 
Lincoln,  Nebraska 
Lubbock,  Texas 
Rapid  City,  South  Dakota 
Tempe,  Arizona 


'Station  Headquarters:  240  W.  Prospect  St.,  Fort  Collins,  CO  80526 


